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Abstract 



We improve our previous QCD sum rule calculation on gKNA and gKNT, 
coupling constants by including the contributions from higher dimensional 
condensates, (qg s a ■ Gq) and (qq)(^-G 2 ), in the OPE. It is found that the 
contribution of these condensates is non-negligible compared to that of the 
quark condensates. Using a best-fit analysis we find \gKNA\ = 2.49 ± 1.25 
and \g K NT,\ = 0.395 ± 0.377. 
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I. INTRODUCTION 



To understand kaon-nuclear physics, it is important to know the hadronic coupling 
strengths involving the kaons. Among them, Qkna and gxNZ are the most relevant cou- 
pling constants. In contrast to g n NN, however, the determination of these kaon couplings 
has some difficulties both in the experimental side and in the theoretical side, e.g. see 0. 

Among other theoretical approaches, QCD sum rule method @-[|] has been used to 
extract these kaon couplings. However, compared to the large number of works devoted to 
9nNN, there have been only few QCD sum rule estimates on gKNA and Qknt, H~@]> f° r which 
there are still ambiguities in among the calculations. Thus the results are quite different 
from each other. More detailed analyses are needed both experimentally and theoretically 
to understand this discrepancy, and to understand kaon-nuclear physics. 

In Ref. the OPE was calculated only up to the leading term coming from the 

quark condensate and to leading order in m s in the sum rule structure proportional to $75. 
However, the next leading term, dimension 5 (qg s a-Gq) may contribute to the OPE side with 
considerable amount as in nucleon mass sum rule |l(J. In addition, operators of dimension 7 
may also be important in the OPE side as a further power correction. Thus, in this paper we 
re-analyze our QCD sum rule calculation including higher dimensional condensates, such as 
(qg s cr ■ Gq) and (qq)(^-G 2 ) , and study the contribution of these condensates on the previous 
results. 

In Sec. H we present our sum rules for gKNA and gKNH, and Sec. |T| we discuss some 
uncertainties in our sum rules and summarize our results. 



II. QCD SUM RULES FOR g KNA AND g KNS 



We will closely follow the procedures given in Refs. JTl],[|||[7| . Consider the three point 
function constructed of the two baryon currents T)b, Vb' and the pseudoscalar meson current 
3s- 

A(p,p',q) = J dxdy{Q\T{ VBI {x) 3 Mn B ^W)e i{p '- x - q - y) . (1) 
In order to obtain gKNA, we will use the following currents for the nucleon and the A |I2|,3- 

Vn = e abc (u^Cj IJ ,Ub)j 5 j' J 'd c , (2) 



Va 



^abc 



{ulC'yuSbhs'fdc - (c^C'7 At s 6 )7 5 7'H 



(3) 



where u and d are the up and down quark fields (a, b and c are color indices), T denotes the 
transpose in Dirac space, and C is the charge conjugation matrix. For the K~ we choose 
the current 



3k- = s^7 5 ^^. 

The general expression for A(p,p', q) has the following form 



(4) 
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A(p,p', q) = F^p^p' 2 , g 2 )i 75 + F 2 (p 2 ,p' 2 , q 2 )<$i~f 5 

+ F 3 (p 2 ,p' 2 , q 2 )Pi l5 + F,(p 2 ,p' 2 , q 2 )^ l5 q,p' u , (5) 



where q = p' —p and P = Recently, in Ref. |T3| it was reported that in the case of g n 



2 



tNN 



the cr /i ^75 structure is independent of the effective models employed in the phenomenological 
side and further provides the 7rNN coupling with less uncertainties from QCD parameters. 
Motivated by this result gKNA and Qknji were calculated from this structure in Refs. 0|§. 
In this paper, however, we construct the sum rule for only the ^75 structure as before, and 
compare this with our previous one. 

On the phenomenological side, keeping the first two terms we have 

x x Mb v 1 f K m 2 K 



(p 2 - M 2 N ){p> 2 - Ml) va "»*™" q 2 - m 2 K 2m, 
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M' B ... > 1 f K m 



+ XnXa * (p* - M 2 N )(p" - Ml)^ l5)9KNA * g 2 - m\ 2m q 
+ higher resonances, (6) 

where M B = \{M N + M A ), and M' B = \(M N - Ml). Here A* means the A (1405), and we 
introduce (-) sign for the A (1405) mass because it is a negative parity state. However, this 
is not relevant in the following calculation. A at, Aa and Aa* are the coupling strengths of 
the baryons to their currents. m q is the average of the quark masses, fx the kaon decay 
constant and m K the kaon mass. We take fx = 0.160 GeV and m s = 0.150 GeV. 

As for the OPE side, the new contribution from the quark-gluon condensates is given by 

- yf^i H-P 2 )((<19 S <7 ■ Gq) + (sg s a ■ Gs)), (7) 
and from dimension 7 ops. 

where we take the limit p' 2 — > p 2 and let (uu) = (dd) = (qq), (ug s a ■ Gu) = (dg s a ■ 
Gd) = (qg s a ■ Gq). Here we collect only the terms which contribute to the <j[/q 2 structure 
such as Figs. [I] and ||. Using the standard values for (sg s a ■ Gs) = 0.8 (qg s o ■ Gq) and 
(qg s (y ■ Gq) = m\ (qq) = 0.8 (qq) |R1 the sum rule after Borel transformation to p 2 = p' 2 
becomes 

X ^ J^fT (e- M » /M2 ~ z~ Ml/M2 ) 9knJ^ + A U~*sl* - e~^\ = 
MX — Mfj v ' 2m q v ' 

2 ( 33 „ , . ,11m 2 21 , „ , , 9 ,m s , . l.a.w,. ,\ , . 

F\M + (— -f - ——)E M 2 + (-i(sa) + -(—G )) (qq). (9) 



V 3 V40vr 2 v 60tt 2 lOOvr 2 ' ' v 3 ' 8 s tt 

Here, A is the unknown constant coming from A a* • gKNA*, and 
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where So is a continuum threshold. One should be cautious, however, that there may be 
non-accounted terms, which can not be neglected by using this simple Borel transformation 



15,16 



N and A 12. 



For \ N and Aa, we use the values obtained from the following baryon sum rules for the 

<M 6 + bE^M 2 + -a 2 = 2{2Txf\ 2 N e- M2 ^ M \ (11) 
3 



£ 2 A M 6 + -am s (l - 3 7 )£ A M 2 + bE^M 2 + -a 2 (3 + 4 7 ) = 2(2vr) 4 A 2 eT M H M \ (12) 
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where a = — (2n) 2 (qq), b = ir 2 {^G 2 ), and 7 = (ss)/(qq) — 1 ~ — 0.2. We use different 
thresholds for X N and Aa in Eqs. ( |TTD and (]T2]). We take = (1.440 GeV) 2 for the nucleon 
sum rule and sa = (1.405 GeV) 2 for the A sum rule considering the next excited nucleon 
and A state, respectively. 

9kna > however, does not display a plateau as a function of the Borel mass. This is 
because there is no usual power correction term like (^2, -pr, and so on) in the r.h.s. of Eq. 
(^) even including up to dimension 7 operators. We need more higher dimensional operators 
to get those terms. Thus, in this case we prefer to use a best-fit method. Eq. (H) has the 
following form : 

9kna ■ fi(M 2 ) + A ■ f 2 (M 2 ) = f 3 (M 2 ). (13) 

Then, we get gxNA and the unknown constant A by minimizing (gxNA ■ fi + A ■ f 2 — fz) 2 
with a fixed so and an appropriate Borel interval: 

/ (9kna ■ fi + A ■ f 2 - / 3 ) dM 2 = minimum. (14) 

Jm 2 - 

mm 

We fix the continuum threshold Sq = 2.074 GeV 2 taking into account the next term from 
the N(1440), i.e., A^(1440) — > A, in the phenomenological side. 

The Borel interval M 2 is restricted by the following conditions : OPE convergence and 
pole dominance. The lower limit of M 2 , M^ in is determined as the value at which the 
contribution of the highest dimensional operators is less than 10% of total OPE. The upper 
limit M^ oa , is determined by restricting the continuum contribution to be less than 50%. 
Then, we get 

\9kna\ = 2.49, 

\A\ = 0.00174 GeV 7 , (15) 

and the Borel interval (0.478, 1.068) GeV 2 for basic inputs (i.e., (qq) = - (0.230 GeV) 3 , 
(^G 2 ) = 0.012 GeV 4 , m s = 0.150 GeV, and m 2 = 0.8 GeV 2 ). Here we denote 
the absolute value because we can not determine signs of the coupling strengths (Ajy, 
Aa and Aa*) in the baryon sum rules. We also calculate the average deviation 5 = 
Ef |1 - RHS(M 2 )/LHS(M 2 )\ /N = 8.8 x 10~ 2 to test the reliability of our fitting, and it 
shows that the deviation is less than 10 %. 
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Table [I] shows variations of Qknk for other input parameters, which are coming from the 
uncertainty of the basic inputs. For example, the first line in Table | shows that \gKN\\ — 
3.11 (or 1.93) if we change the quark condensate to (qq) = -(0.210 GeV) 3 (or -(0.250 GeV) 3 ) 
while other basic inputs are fixed. In the last line we take j«q = 0.6 GeV 2 from the lowest 
value of the standard QCD sum rule estimate |1~4[], and 1.4 GeV 2 which was evaluated in 
the instanton vacuum in Ref. |I7|. Total variation is about ± 1.25 on the above cjkna value. 



On the other hand, the unknown constant \A\ varies from 0.00120 to 0.00203 GeV 7 . 

Next, consider cjkny,- The current of S° is obtained by making an SU(3) rotation from 
the nucleon current 



^ = ^2 e ahc [(^(77^)757^4 + ((^C^s b ) 757 M «c] • (16) 
In this case the contribution of the quark-gluon condensate is given by 

- V2^\n(-p 2 )((qg s a ■ Gq) + {Sg.tr ■ Gs}), (17) 
and from dimension 7 ops. 

+ v / 2^^((^) + (-))(fG 2 ). (18) 
Then, within the same approximation as before we get the following sum rule. 

Xn ^m 2 -m 2 i e ~ e ) 9KN *^T + B \ e ~ p - 



where Mb = \{Mn + Ms) and N* is N(1440). B is the unknown constant coming from 
Aat* ■ gxN*Y,- Again for As, we take the value from the following sum rule for the £ flT2"|,|3"|l : 

EfM 6 - 2am s (l + -f)E*M 2 + bE^M 2 + -a 2 = 2(27r) 4 A|e- M ^ M2 . (20) 

3 

We fix the continuum threshold ss = (1.660 GeV) 2 considering the next £ state, £ (1660). 

Using the continuum threshold Sq = 2.356 GeV 2 taking into account the next term from 
the N(1535), i.e., A^(1535) — > £, in the phenomenological side we get 

\9kns\ = 0.395, 

\B\ = 0.00148 GeV 7 (21) 

for the same basic inputs. The Borel interval is (0.488, 1.584) GeV 2 and the average deviation 
of the fit 5 is 9.7 % in this case. We present the variation of gKNY, on other parameters in 
Table |Il|. The total variation is about ± 0.377. On the other hand, \B\ varies from 0.00117 
to 0.00184 GeV 7 . 
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III. DISCUSSION 



SU(3) symmetry, using de Swart's convention predicts 



1 

gKNA — — — 2a D )g 7rNN , 

9knt, = + (2«d - l)g nN N, (22) 



where old is the fraction of the D type coupling, a, D — S~f- m Table [TTJ we compare our 



results with previous QCD sum rule estimates |||8]|| and an SU(3) symmetry prediction, 
where we denote the error-bar allowing for SU(3) symmetry breaking at the 20 % level. 
Here we take from a recent analysis of hyperon semi-leptonic decay data by Ratcliffe, 



«£>=(]. 64 [pQ| ; and g^NN from an analysis of the np data by Ericson et al. pl| , p^jvat = 13.43. 
A comparison to fitting analyses of experimental data |E| is also provided. SU(3) symmetry 
predicts | gKNA / gKNY. \ = 3.55 taking ao = 0.64, while our results show that this ratio is 6.30 
using the basic inputs, and the order of SU(3) symmetry breaking is rather huge. 

Let us remark on g n NN which was calculated in Ref. flTT|j3| using the three-point function 



method. After including dimension 5 and 7 condensates as in the previous section the sum 
rule becomes 



~ (i^ 4 " ^oM 2 + ^G 2 )) (ft), (23) 

where C is the unknown constant from Atv» • g n NN* and f n = 0.133 GeV. The contribution 
of the quark-gluon condensates in the OPE side is important as in the gxNA and gxNY, 
sum rules. In this case we use the PCAC relation f^m^ = —Am q (qq) first, then the quark 
condensate becomes an overall factor on both sides. However, the coupling strength is 
still related to the quark condensate as shown in Eq. ([TT]) . 

Following the same method in the previous section, and using (qq) = -(0.230 GeV) 3 , 
(^-G 2 ) = 0.012 GeV 4 , and s = 2.074 GeV 2 as a pure continuum threshold we get 

|#7rAw| = 3.65 ± 2.31, 

\C\ = 0.00261 ± 0.00091 GeV 7 , (24) 

the Borel interval (0.460, 1.110) GeV 2 , and the average deviation of the fit 5 9.3 % at the 
central value. Here the uncertainty comes from using different input parameters, i.e. (qq) — 
-(0.210 GeV) 3 (or -(0.250 GeV) 3 ) , (^G 2 ) = 0.015 GeV 4 , and m 2 = 0.6 (or 1.4 GeV 2 ) as 
before. In this case the most error bar comes from uncertainties of the quark condensate, 
i.e. from the coupling strength A^. 

Now, let us discuss some uncertainties in our sum rules. In Eqs. (|^), (f03) and (E3) the 
contribution of the quark-gluon condensate is about 25 %, 40 %, and 20 %, respectively, of 
the leading term at M 2 = 1 GeV 2 . Thus the accurate value of this condensate is one of 
important factors in our sum rules, and a more precise estimate may be needed (e.g., see 
Ref. ra. 
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As we mentioned before, we need more higher dimensional operators to get some power 
correction terms in our sum rules. Their contribution will be much smaller than that of 
dimension 7 operators at the relevant Borel region around M 2 ~ 1 GeV 2 . However, those 
operators may contribute because the lower limit of the Borel interval for each coupling 
constant is much less than 1 GeV 2 in our sum rules. 

We find that the coupling constants become 2 or 3 times larger than the previous ones if 
we take the coupling strengths (Ajy, Aa and As) from the chiral-odd baryon sum rules |]T2"|,|3"|] . 
For example, we get 7.04, 0.890, and 14.49 for \qkna\, \9kns\, an d \9ttNn\i respectively, for 
the basic inputs. Because the coupling strengths from each baryon sum rule (the chiral-even 
and chiral-odd) are not the same in the whole Borel region and the discrepancy between 
the coupling strengths is larger in the low Borel region, we get quite different coupling 
constants. Of course, it should be judged by the stability of the sum rule whether one 
chooses the coupling strengths from the chiral-even sum rules or those from the chiral-odd 
sum rules. 

As a final remark, in the case of g^NN it was shown that there is a higher pseudoscalar 
resonance contamination from the tt (1300) and it (1800) in the three-point function method 
p4| . Maybe there is a similar contamination from the K(1460) and K(1830) |25j on the kaon- 



baryon couplings. Although the masses of the K(1460) and K(1830) are quite uncertain and 
these states need further experimental confirmation, we can briefly estimate the contribution 
of the K(1460) as done in Ref. |f24| . Using the parameters from recent works p6 , we get 



Q 2 + m 2 



M 



= 21.3 and 2.2 MeV (25) 

=1 GeV 2 



for the kaon and K(1460), respectively. Here Jm is the decay constant and tum is the meson 
mass. We take fx = 108 MeV, fxium) — 3.3 MeV and m K = 496 MeV, rriK{ims) — 

1.45 

GeV in Ref. [ 26f . Comparing the values in Eq. (|25|) to those for the pion and ir (1300) ||26|| . 



i.e. 1.7 and 0.4 MeV, the contamination from the excited kaon state on the kaon-baryon 
couplings seems smaller than that from the excited pion state on q-knn. 

In summary, including higher dimensional condensates we re- analyze our previous QCD 
sum rule estimate on Qknk and gxNT, in the ^75 structure. The contribution of dimension 
5 quark-gluon condensates is comparable to that of the leading term, and the present result 
is much different from the previous one. 
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TABLES 

TABLE I. gxNA and its variations. Other inputs mean other possible inputs coming from the 
uncertainty of the basic inputs. 



basic inputs 


other inputs 


variations 


(qq) = -(0.230 GeV) 3 


-(0.210 GeV) 3 , -(0.250GeV) 3 


+ 0.62, - 0.53 


\-^-(_x ) = U.Uiz Lxev 


U.UlO LxeV 


+ U.40 


m s — u.iou Lxev 


n 1 on n 1 cn r^^AT 
U.lzU, U.loU Lxev 


n Q£ i n Qfi 
— U.oO, + U.oU 


Trig = 0.8 Lxev 


O.o, 1.4 Lrev 


- u.iy, o.oz 


TABLE II. 


9k NT, and its variations. The same as in Table |. 




basic inputs 


other inputs 


variations 


(qq) = -(0.230 GeV) 3 


-(0.210 GeV) 3 , -(0.250 GeV) 3 


+ 0.057, - 0.061 


(^G 2 ) = 0.012 GeV 4 


0.015 GeV 4 


+ 0.201 


m s = 0.150 GeV 


0.120, 0.180 GeV 


- 0.084, + 0.086 


ml = 0.8 GeV 2 


0.6, 1.4 GeV 2 


+ 0.067, - 0.198 



TABLE III. Comparison of coupling constants. 



Sources 



9KNA 



9KNY, 



SU(3) with 20 % breaking 



Experimental fitting [22] 
Ref. § 
Ref. 1 
Ref. @ 
Present work 



16.0 ~ - 10.7 

- 13.7 

10 ± 6 
2.37 ± 0.09 

10 ± 2 
2.49 ± 1.25 



3.0 ~ 4.5 

3.9 
3.6 ± 2 
0.025 ± 0.015 

0.75 ± 0.15 
0.395 ± 0.377 
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FIGURES 



FIG. 1. Contribution of dimension 5 operators. The solid lines are quark propagators and the 
wavy line is a gluon propagator. The dotted line denotes a meson. 

FIG. 2. Contribution of dimension 7 operators. The same as in Fig. |l| 
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